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Abstract 

We consider the open spin-s XXZ quantum spin chain with nondiagonal boundary 
terms. By exploiting certain functional relations at roots of unity, we derive a gen- 
eralized form of T-Q relation involving more than one independent Q(u), which we 
use to propose the Bethe-ansatz-type expressions for the eigenvalues of the transfer 
matrix. At most two of the boundary parameters are set to be arbitrary and the bulk 
anisotropy parameter has values rj = ^ , ^ , . . . . We also provide numerical evidence 
for the completeness of the Bethe-ansatz-type solutions derived, using s = 1 case as an 
example. 



1 cmail:ribaiyas@svsu.edu 
2 e-mail: rmurgan@svsu.edu 



1 Introduction 



Considerable effort has been put into solving integrable quantum spin chains for many years. 
In particular, integrable open quantum spin chains have attracted much interest over the 
years. In this regard, open XXX and XXZ quantum spin chains have been extensively 
investigated due to their growing applications in fields of physics such as statistical mechanics, 
string theory and condensed matter physics. Much progress has been made on the topic up 
to this point. Numerous successes in the past [l]-[8] (also refer to [9]- [31] and references 
therein, for other related work on the subject) have motivated further investigations of these 
models. In addition, in a series of publication, Bethe ansatz solutions have been derived 
for open spin- 1/2 XXZ quantum spin chain where the boundary parameters obey a certain 
constraint. Readers are urged to refer to [32]-j39] for related work on the subject. Two sets of 
Bethe ansatz equations are needed there to obtain all 2 N eigenvalues, where N is the number 
of sites. A special case of the above solution was generalized to open XXZ quantum spin 
chain with alternating spins by Doikou [10] using the functional relation approach proposed 
by Nepomechie in [31] to solve the spin-1/2 case. In [H], related work was carried out using 
the method in [32]. In [12], the spin-1/2 XXZ Bethe ansatz solution (for boundary parameters 
obeying certain constraint) is generalized to the spin-s case by utilizing an approach based 
on the Q-operator and the T-Q equation [43] (see below), which was developed earlier for the 
spin-1/2 XXZ chain in [38] and subsequently applied to the spin-1/2 XYZ chain in [2]. Two 
sets of Bethe ansatz equations are also needed there to produce all (2s + 1)^ eigenvalues, 
where again N represents the number of sites. This was later followed by another work for 
spin-s with such constraint removed, but limiting the Bethe ansatz solutions for cases with 
at most two arbitrary boundary parameters for some special values of the bulk anisotropy 
parameter [35], namely r\ = ^y, with p being even integers |45j . 

In a number of works cited above, the well known Baxter T — Q relation [13], with the 
following schematic form 

t(u)Q(u) = Q(v) + Q{w) (1.1) 

has provided a way to obtain the Bethe ansatz equations for the eigenvalues of the transfer 
matrix t(u). In [4~6| Wf\ . a generalization of this relation involving two Q(u) of the following 
form for the open spin-1/2 XXZ quantum spin chain was given: 

t{u)Q x {u) = Q 2 (v) + Q 2 (w) 

t(u)Q 2 (u)=Q 1 (v') + Q 1 (w') (1.2) 

Motivated by this solution, in this paper, we obtain the corresponding solution for the 
open spin-s XXZ quantum spin chain. In addition, our work is motivated by the fact that 
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such T — Q relations for an open spin-s XXZ quantum spin chain are novel structures and 
therefore merit further studies and investigation. We remark that a more general form of 
T—Q relations were found in [18] involving multiple Q(u) functions. Moreover, the relation of 
s — 1 case to the supersymmetric sine-Gordon (SSG) model [H]-[5l], especially the boundary 
SSG model [55]- [59], has inspired us to consider the problem. We stress that these results 
hold for cases with at most two arbitrary boundary parameters at roots of unity, namely 

when the bulk anisotropy parameter has vales rj = y , ^ , We follow similar approach 

as given in [31]- [37] and [16] that was used to solve the s = 1/2 case, which is based on 
functional relations obeyed by transfer matrix at roots of unity. This yields Bethe-ansatz- 
type solutions which give the eigenvalues. Our numerical analysis for s — 1 case for N = 2 
and p = 3,5 yields all the eigenvalues as given in Tables 1 and 2. As in |42j, we rely on 
fusion [5], [60]- [66], the truncation of the fusion hierarchy at roots of unity [67]- [69] and the 
Bazhanov-Reshetikin [70] solution of the RSOS models. 

The outline of the paper is as follows: In Sec. 2, we review the construction of the fused 
R [60]-[64], [7_l]-[75] and K T [5], [651 [66] matrices from the corresponding spin-1/2 matrices. 
One can refer to [761 E7J for some original work on spin-1/2 K T matrices. We then review 
the construction of commuting transfer matrices from these fused matrices (using Sklyanin's 
work [4], which relies on Cherednik's previous results |80j), together with some of their 
properties. Fusion hierachy and functional relations obeyed by transfer matrices are also 
reviewed. In Sec. 3, the generalized T — Q relations are given along with some arguments 
behind their structure. This is done by exploiting the reviewed functional relations obeyed 
by the transfer matrices. From this, we derive the Bethe-ansatz-type equations for cases 

with at most two arbitrary boundary parameters at roots of unity, e.g. r\ = y , ^ , 

We then present numerical results in Sec. 4 to illustrate the completeness of our solution, 
using s = 1 as an example. Here, the Bethe roots and energy eigenvalues derived from the 
Bethe-ansatz-type equations (for some values of p and N) are given. We remark that these 
eigenvalues coincide with the ones obtained from direct diagonalization of the open spin-1 
XXZ chain Hamiltonian. Finally, we conclude the paper with discussion of the results and 
potential future works in Sec. 5. 



2 Commuting spin- 5 transfer matrices and functional 
relations at roots of unity 

In this section, in order to make the paper relatively self contained, we review some crucial 
concepts on the construction of commuting transfer matrices for iV-site open spin-s XXZ 
quantum spin chain. Materials reviewed here on fused R, K T and higher spin transfer 
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matrices are mainly reproduced from [32] • Like the commuting transfer matrix for s = 
1/2, constructed in [3], which we denote (following notations adopted in [32]) by t^'^(u), 
whose auxiliary space as well as each of its N quantum spaces are two-dimensional, one can 
construct a transfer matrix $ ,s >{u) whose auxiliary space is spin-j ((2j' + l)-dimensional) and 
each of its N quantum spaces are spin-s ((2s + l)-dimensional), for any j,s G {|, 1, |, . . .} 
using the fused R [60]- [63], [H]-[75] and [5], (65J [66] matrices. These matrices serve as 
building blocks in the construction of the commuting transfer matrices for higher spins. We 
list them below with some of their properties. The fused- R matrices can be constructed as 
given below, 



2j 2s 



p w p w n n r{ z£ ] ( u +(k+i-j-s-i) V ) 



{b} 



(2.1) 



k=l 1=1 



where {a} = {a 1; . . . , a 2 j}, {b} = {bi, . . . , &2s} ; and Pf-, is the symmetric projector given by 



1 2J ( k 
P {a} = (2j)l II I J2 V < 



(2.2) 



k=i \i=i 

V is the permutation operator, with V ak ,a k = 1; similar definition also holds for P i 



R^^)[u) is given by 

R^'^(u) 



{by 



\ 



(2.3) 



/ sh(w + rj) 

shw sh rj 

sh?7 shw 
sh{u + r]) ) 

where rj is the bulk anisotropy parameter. The fused R matrices satisfy the Yang-Baxter 
equations [78j [79] 



V 



Ri% } ^-v)R^ {c} (u)R^% } (v) 



>(i,s) 



(2.4) 



The construction of the fused K matrices now readily follows [5], [651 66J, 



M 



2j 

n 

fc=i 



fc-i 



1=1 



x /^( u + (fc-j-_), 7 )jp+ } , 

where -^"^(w) is the 2x2 matrix whose components are given by [761 177] 
-^Ti ( M ) = 2 (sh a_ ch /3_ ch u + ch a_ sh /3_ sh ii) 
(u) = 2 (sh «_ ch /3_ ch -u — ch a_ sh /3_ sh u) 
K{ 2 (u) = e e ~ sh 2u , #2i(«) = e ~ 9 ~ sh 2u , 



(2.5) 



(2.6) 
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where «_,/?_, are the boundary parameters. The fused K matrices satisfy the boundary 
Yang-Baxter equations [HQ] 



{6} V / {a}{6}^ ' {a} V I {a}{b} 

In addition, the fused K + matrices are given by 



,CmO 



-0)/ 



>0» 



M — V) 



with the following normalization factor 



(a_,/?_,0_)-K-a + ,-/3 + ,0+) 



2J-1 / 

= n n-^ + c + k + 1 - 2 ^»] • 
i=i fc=i 



From the fused matrices, one constructs the higher spin transfer matrix t^' s \ 

t^\u) = tr {a} K$\u) T?*\u) K-^(u) f^\u) . 
The monodromy matrices are given by products of N fused R matrices, 

(«), 



u 



(2.7) 



(2i 



(2.9) 



(2.10) 



p(i> s ) /„, \ p0'. s ) 



>0» , 
l {a},{feW} ( 



>0» / 
l {a},{bl N ]} [ 



(2.11) 



These transfer matrices commute for different values of spectral parameter for any j,j' G 
{§, 1, |, . . .} and any s G {|, !,§,.. .}, 







(2.12) 



Furthermore, they also obey the fusion hierarchy [5j H2J [651 66]j 



_ t(i«)( M ) = *(?.«) ( u _ (j _ 1)^) + _ jy) - (j + i)^) > (2-13) 



j = 1, |, . . ., where t^°' s ^ = 1, and S^(u) is given by 



u 



(2.14) 



where 



2s-l ^ 

f| e(« + (s - fe + ^) 



n 2JV 



fc=0 



sh(2w) sh(2« + 477) 



sh(2w + 77) sh(2u + 3t?) 

(u) = 2 4 sh(w + a_ + 77) sh(w — a_ + 77) ch(-u + /3_ + 77) ch(w — /3_ + 77) 

x sh(u + a + + 77) sh(?j — ct + + 77) ch(u + /3+ + 77) ch(u — (3 + + rj) . (2.15) 



1 See the appendix in [42] for more details on the fusion hierachy. 



To avoid confusion, we emphasize that the 5^ (u) in [42j differs from the one given here by 
a shift in 77. 

Next, we list a few important properties of the rescaled "fundamental" transfer matrix 
t^' s \u) (defined below), which are useful in determining its eigenvalues. Following the 
definition of t^ ,s \u) as in [12], we have 

= (1 * * ( *' a) (")> (2- 16 ) 

g\2' s )[uY N 

where 

2s- 1 j 

= JJ s h(w + (s - fc + -)//) . (2.17) 

k=l 

This transfer matrix has the following useful properties: 

fth s \ u + ztt) = t(i s )( M ) (ztt - periodicity) (2.18) 

- n) = t^' s \u) (crossing) (2.19) 
t(i s )(0) = -2 3 sh 2iV ((s + ^)r/)ch?7sha_ch/3_sha + ch/3 + I (initial condition) (2.20) 

t (1 2' s \u) = 2 3 sh 27V M[-sha_ch/3_sha+ch/3+ch 2 M 

T)=0 

+ ch a_ sh /3_ ch a + sh /3 + sh 2 u 

- ch(fl_ - # + )sh 2 uch 2 w I (semi-classical limit) (2.21) 

where I is the identity matrix. 

Due to the commutativity property (12.121) . the corresponding simultaneous eigenvectors 
are independent of the spectral parameter. Hence, ( 12.181) - ( 12.211) hold for the corresponding 
eigenvalues as well. In addition to the above mentioned properties, for bulk anisotropy 
parameter values r] = ^y, with p = 1,2,..., the "fundamental" transfer matrix, t^' s \u) 
(and hence each of the corresponding eigenvalues, A^' s \u)) obeys functional relations of 
order p + 1 [34|-[36j, 

t(M( M )t(^)( M + v ) . . . t^' s \u + prj) 

- 5 {s \u-r])t^' s \u + r])t^' s) {u + 2 V )...t^' s) {u + {p-l)r]) 

- 5 {s \u)t^' s) {u + 2ri)t^' s \u + 377) . . .t^' s \u + prj) 
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- 5 {s \u + rj)t^' s \u)t^' s) (u + 3r])t^' s) (u + 4rj)... t^' s) (u + prj) 

- 5 {s \u + 2r])t^ s \u)t^' s \u + r])t^' s \u + At]) . . . t&'\u + prj) - . . . 

- 5 (s \u + (p - l)ri)t^' s \u)t^' s \u + 77) . . . t^' s \u + (p - 2)77) 

+ ... = /(u). (2.22) 

The scalar function /(tt) (which can be expressed as /(tt) = fo(u)fi(u)) is given in terms of 
the boundary parameters a T , f3 T , 9 T (for odd p) by 



f (_l)*+l 2 -4*iV sh 4s7V ((p + th 2 ((p + 1)«) 

„ _ 1 3 5 



/oH = < 



2 ' 2 ' 2 ' • • • 



.^TV+^-^piV ch 4s7V ^ + th 2 ^ + ^ 



u 



(2.23) 



S = l,2,3,... 



and 

/i(«) = "2 3 - 2p ( 

ch ((p + l)a_) ch ((p + ch ((p + l)a+) ch ((p + l)/3+) sh 2 ((p + 

- sh ((p + l)a_) sh ((p + sh ((p + l)a+) sh ((p + l)/3+) ch 2 ((p + l)u) 

+ (-l) Jv ch((p + l)(^-^+))sh 2 ((p + l)u)ch 2 ((p + l)n)) , (2.24) 

for s = I , § , I . . . and 

fi(u) = (-l) JV + 1 2 s -*( 

ch ((p + l)a_) ch ((p + ch ((p + l)a+) ch ((p + l)/3+) sh 2 ((p + l)u) 

- sh ((p + l)a_) sh ((p + sh ((p + l)a+) sh ((p + l)/3+) ch 2 ((p + 

+ ch ((p + 1) (#_-#+)) sh 2 ((p + 1)m) ch 2 ((p + 1)m) ) , (2.25) 

for s = 1 , 2 , 3 . . .. Note that f(u) satisfies 

f(u + rj) = f(u), f(-u) = f(u), (2.26) 



and 



fo(u) 2 = f[S { s \u + 3V ) 



(2.27) 



3=0 



where Sq\u) is given by ( 12. 15ft . 



3 Generalized T-Q relations and Bethe ansatz 



In this section, we give the main results of this paper. We derive the generalized T — Q 
relations for the transfer matrix eigenvalues and obtain the Bethe-ansatz-type equations, for 
cases where at most two of the boundary parameters a± or (3± are arbitrary, by adopting 
the steps given in [46] while setting 9_ = 9 + = 9, where 9 is also arbitrary. More on this is 
given below. 



3.1 T — Q relations 

The transfer matrix t^' s \u) and its eigenvalues (A^' s \u)) obey the functional relations 
( I2.22j) . We exploit this fact to obtain the T — Q relations. Following [70], one could recast 
the functional relations as the condition that the determinant of a certain matrix vanishes, 
namely 



det M(u) = 0, 
where Ai(u) is given by the (p + 1) x (p + 1) matrix 



M(u) 













(3.1) 







\ -/i (2) (m - 77) 



-h^{u + (p- 1)77) A^' s \u + pr]) J 



(3-2) 



where (u) and (u) are functions which are z7r-periodic, but otherwise not yet specified. 
We note that the above matrix has the following symmetry, 



TM(u)T- 1 = M(u + 2r]), T = S 2 



(3.3) 



where S is given by, 



/ 1 .. 
1.. 



S 



0.. 

\ 1 .. 



\ 




1 
0/ 



(3.4) 



Assuming that 



det M(u) = 



(3.5) 
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then M.(u) has a null eigenvector v(u), 

M(u)v(u) = 0. (3.6) 
The symmetry (I3.3P is consistent with 

T v(u) = v(u + 2rf) , (3.7) 
which implies that v (u) has the form 

v(u) = (Q^u) , Q 2 {u + 77) , . . . , Q 1 (u - 2rj) , Q 2 {u - r,)) , (3.8) 

with 

Qt(u) = Q^u + in) , Q 2 {u) = Q 2 {u + m) . (3.9) 
That is, the components of v (u) are determined by two independent functions, Qi(u) and 

The null eigenvector condition ( 13. 6 j) together with the explicit forms of A^(w) and 
given by (13. 2 p and (13. 8p respectively, now lead to the following T — Q relations, 

5(°\u)Q 2 (u + V ) 5^(u- V )Q 2 (u- V ) 
{) hW(u) Q x {u) ^hMiu-ri) Q 1 (u) ' 1 } 

= h^u-v)^^ + h^u)9^. (3.11) 

Due to the crossing symmetry fl 2 . 1 9 [) and 

5^{u) = 5 {s) {-u-2r]), (3.12) 

which is the crossing property for 5^ s \u), it is then natural to have the two terms in (13.101) 
transform into each other under crossing. Hence, we set 

hP-)(u) = h,W(-u -2rj), (3.13) 

and we make the following ansatz 

Mi 

Qj(u) = Y[sinh(u - u^ } ) sinh(u + +77), (3.14) 
fc=i 

which is consistent with the required periodicity (13. 9p and crossing properties 

Qj(u) = Qji-u-rj), (3.15) 



2 In [45], all of the matrices M.(u) possess a stronger symmetry, S M{u)S 1 = A4(u + 77), implying the 
null eigenvector with single Q(u). We refer the reader to Sec. 3 of [46] for more detail discussion on this. 



8 



where j = 1,2. In (I3.14p . {u^} represents the Bethe roots (or zeros of Qj{u)) and there are 
Mj of these roots. Further, one can verify that the condition detj^t(u) = indeed implies 
the functional relations (12.221) . if w(u) satisfies 



p-i p 

f{u) = w{u) \\ tfW(« + ^) + _ ] | 6M(u + jrj), 
- LJ - w(u) - LJ - 



J'=0,2 V 



(3.16) 



i=i,3, 



where 



Will) = 



nj ,,,.../' 2, (" • ./'/) 



(3.17) 



It follows from (I3.16P that the process of finding w(u) reduces to solving a quadratic equation, 
which when used together with (I3.13P and (I3.17p . yields the explicit form for the function 
h^(u). Here, we consider even number of sites, N. Below, we give the solutions of (I3.17P 
for h^(u) for two cases: 

I. /?_ and (3 + arbitrary while setting a± = 0, 6L = 9 + = 9 = arbitrary. 



h {1 \u) = 4 



2s-l 



Y[ sh(u + (s - k + -)r]) 



k=0 



2N 



sh 2 (ii + rf) sh(2« + Ai]) 
sh{2u + 3rf) 



Mi = sN + -(p + l), M 2 = M 1 -1. 



(3.18) 



II. a_ and a + arbitrary while setting (3± = 0, 9_ = 9 + = 9 = arbitrary. 



■it 



2s-l 



-I 2N 



J J sh(u + (s - k + -)rj) 



k=0 



ch 2 (u + rj) sh(2ii + Arf) 
sh(2w + 3r]) 



M 1 = sN+~(p + l), M 2 = Mx-l. 



(3.19) 



Now, using the analyticity of A^2' s \u), given by (I3.10P and (13. lip , one can write down 



the Bethe-ansatz-type equations for the zeros {uf ,uf} of Q\{u) and Q2{u) 



S^uP) hV>(uf>-7i) 



(1) 



h^(uf -rj) 



(2)> 



Q2(uf ) - V ) 
Q2(uf ) + V ) 

Qi(uf+v) 
Qi{uf-rj) 



j = l,2 } .. . ,M a 
j = l,2,... ,M 2 



(3.20) 
(3.21) 



We remark here that for each case, there are more than one solutions for h^ 1 ' (u) that corre- 
spond to the above expression for w(u). The solutions found are largely by trial and error, 
verifying numerically for small values of N that the eigenvalues can indeed be expressed as 
(EHUD , (I3TTD with Q(u)'s of the form f[3TT4) . 

To summarize, we have proposed that for the case where the bulk anisotropy parameter, 
i] = -^j with p being odd integers and that at most two of the boundary parameters are 
arbitrary, the eigenvalues A^' s \u) of the transfer matrix t^' s '(u) for two cases (I and II) are 
given by a generalized form of T — Q relations ( 13. 101) . (13.111) . with Q±(u) and Q2{u) given by 
(gUD and h {2 \u) given by fl3TT3|) . The h^(u) is given by (ETTgjl and fl3TT9|) respectively, for 
the two cases considered. The zeros {uf ,uf} of Q x {u) and Q 2 (u) are indeed the solutions 
of the Bethe-ansatz-type equations, (13.201) and ( 13.21D . These equations reproduce results in 
[IHJ H7] for s = |. In the following section, we shall use these results (specifically A^2- s )(m) 
which represents the eigenvalues of the rescaled "fundamental" transfer matrix given by 
(I2.16P ) to derive expressions for energy eigenvalues for the case s = 1. 

4 Energy eigenvalues and Bethe roots 

In this section, we provide numerical evidence for the completeness of the Bethe-ansatz-type 
solutions derived in Sec. 3 using the case s = 1 as an example. We derive an expression 
for the energy eigenvalues for the open spin-1 XXZ quantum spin chain and compute the 
complete energy levels for this case from the Bethe roots given by (I3.20p and (I3.2ip . We do 
this for both cases, I and II. 

4.1 Open spin-1 XXZ quantum spin chain: The Hamiltonian 

In this section, we review the integrable Hamiltonian for the open spin-1 XXZ quantum spin 
chain (adopting notations used in |42j). The Hamiltonian is given by 

N-l 

U=Ys H n,n+l + H b , (4.1) 
n=l 

where H n>n+ i represents the bulk terms. Explicitly, these terms are given by [83] . 

H n , n+1 = a n - K) 2 + 2 sh 2 V [< + (S z n ) 2 + (S* n+1 ) 2 - «) 2 ] 

- 4sh 2 (|)(^< + <^) , (4.2) 

where 

o"n = S n ■ S n +i , = S^S^+i + S^S y n+1 , a z n = S^S z n+1 , (4.3) 
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— * 

and S are the su(2) spin-1 generators. Hf, represents the boundary terms with the following 
form (see e.g., [221183] ) 



H b = cuiSl) 2 + a 2 S{ + a 3 {S+) 2 + a 4 {S^) 2 + a 5 S+ S{ + a & S{ 5f 
+ a 7 Sl + a 8 Si Si + (aj o bj and 1 -B- iV) , 



(4.4) 



where S^ = S x ±iS y . The coefficients {a«} of the boundary terms at site 1 are functions of 
the boundary parameters (a_, 6L) and the bulk anisotropy parameter 77. They are given 

by, 



«i 

a i 
a 5 

«7 
«8 



a (ch 2a_ — ch 2/3_ + ch r/) sh 2r/ sh 77 . 
ao sh 2a_ sh 2/3_ sh 2rj , 



aoe 261 " sh 2?7 sh 77 



— -aoe 2e ~ sh 2?7 sh 77 , 
8 

/ 77 77 \ 3 

aoe 9- ^ch /3_ sh a_ ch — + ch a_ sh /3_ sh — J sh 77 ch 5 77 , 
a e~ e ~ I ch /3_ sh a_ ch — + ch a_ sh /3_ sh — J sh 77 ch 5 77 , 



77 



where 



a 



— a e " ( ch /3_ sh a_ ch — — ch a_ sh /?_ sh — J sh 77 ch 2 77 , 

Tj Tj \ 3 

— " ( ch/3_ sha_ ch — — ch a_ sh /3_ sh — J shr/ ch 2 77 , 



sh(a_ - j) sh(«_ + 1) ch(/3_ - |) ch(/3_ + |; 



(4.5) 



(4.6) 



Similarly, the coefficients {b{\ of the boundary terms at site N which are functions of the 
boundary parameters (a+, 0+, 9+) and 77, are given by the following correspondence, 



h = a, 



(4.7) 



The Hamiltonian "H (14 .11) . according to [1], is related to the first derivative of the spin-1 
transfer matrix, namely t^ l,1 \u), which one constructs from t^'^^u) by using the fusion 
hierarchy formula (12.131) . 
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where 5^\u) is given by fl2.14p - p.15p with s = 1. Following [42], we work with the rescaled 
transfer matrix given by 

[shu sh(u + ?7jp v 

where t*- 1 ' 1 ) 9 *(w) is the transfer matrix constructed from "gauge" -transformed R^^u) and 
K T ^\u) matrices [§. We note here that the rescaled transfer matrix does not vanish at u — 0. 

The Hamiltonian "H (14. ip . can now be expressed in terms of the first derivative of 



^ = c«Af(i,i) 9 t( w ) +c (i) I; (4.io) 

GSM 



M=0 



where 

cS 1} = chr/{l6[sh2r7shr/] 27V sh3r/sh(a_ - ^)sh(a_ + ^)ch(/3_ - ^)ch(/3_ + 



x sh(a + -|)sh(a + + |)ch(/3 + -|)ch(/3 + + |)} _i (4.11) 



and 



4 1} = -^&chr?-(iV-l)(4 + ch2r7) + 2;Vch 2 7/ 



sh rj 



2ch2a+^ch?7(3 + 7ch2?7 + ch4?y) + ch2/3 + (4 + 5ch2r? + 2ch4r/)j 



2c/ 

+ 2 clinch 2/3+ (3 + 7ch2?7 + ch4r/) + chr/(5 + 3 ch2r? + 3 ch4r?) j J 
- ch2/3 + (2 + 4ch?7ch377) + chr](5ch2r] + ch4?y) - 2ch2a + (l + ch2?7 

+ ch2/3 + (ch77 + 2ch37/) + ch4r/)j. (4.12) 
In (I4.12p . b and d are given by 

b = 2( - ch2/3_ - ch 3 r/ + ch2a_(l + ch2/3_ ch^)) (4.13) 

and 

d = -4sh3r ? sh(« + + |) sh(a + - |) ch(/3+ + |) ch(/3 + - |) . (4.14) 



3 Such a transformation results in a more symmetric form of these matrices. For a detailed discussion on 
this, refer to Sec. 4 of 142) . 
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4.2 Open spin-1 XXZ quantum spin chain: Energy eigenvalues 



Next, we proceed to the eigenvalues of the Hamiltonian ( I4.10p . Note that f)4.10p implies the 
following result for the corresponding eigenvalues, 

E = c?-^A^ 9t {u) +c W ; (4.15) 
du u=o 

where A^ 1,1 ) 9t {u) represents the transfer matrix eigenvalues which assume the following form 
after using (12TT6|) . (14~8|) and (1479]) . 

- (4-16) 

where 5^(tt) is given by (I2.14p - fl2.15p with s — 1. Furthermore, we have also used the fact 
that A (1,1) = A^\u). Finally, from (EUTJ, fjHD and fO]l . we obtain the energy in 
terms of Bethe roots {m^} for cases I and II. Below, we present the analytic forms of the 
energy eigenvalues for these two cases: 

Case I. (3- and (3 + arbitrary while setting a± = 0, Q- = 9 + = 9 = arbitrary : 

Mi 

E = sh\2 V )J2 ~ ^-^7 ~ (I) +2sh2r / [(iV + l)cthr ? -cth r? ' 
fe=1 sh«' + f ) sh(«^ - f) 



« + ?)ah(««-|) " 2 J 

+ ci iJ C'(0) + 4 ij (4.17) 

Case II. ct_ and ct + arbitrary while setting /3± = 0, 6L = 9 + = 9 = arbitrary : 

Mi 

£ = sh\2 V )J2 ~, (i) 3 ^ w (i) „ +2sh27y[(iV+l)cth7y-th|] 
ttsh(4 ) + f)sh(4 ) -f) 2 

+ ci 1) C(0) + 4 1) (4-18) 
where in f l4TT7|) and fOB"]) . 

C(«) = + % ^(«,- 2), (4.19) 

We recall that Mi = JV + |(p + 1) and M 2 = M l - 1. Note that in (gT7J and (JUSD, 
the contribution to E comes only from {u^} and not from both {uj. } and {u k } as one 
may initially expect. Perhaps, if one uses f l3.10p instead of ( 13. lip in the derivation of E, 
an equivalent expression involving only 2 } or both {i^. } and {u^} may result. This 
however will not affect the numerical value of the energy eigenvalues tabulated in the next 
section. 
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4.3 Numerical results 



In this section, we tabulate the energies computed using (I4.17P and (14. 18j) for some values 
of N ,p (therefore rj) and the boundary parameters {a± , fi± , 9±] with the Bethe roots, 
{u^} } in Tables 1 and 2 for cases I and II respectively. These Bethe roots are obtained 
using McCoy's method [EH [82]. These numerical results demonstrate the completeness of 
the Bethe- ansatz- type equations, ( 13.201) and ( 13.211) for the case s — 1. We checked these 
solutions for chains of length up to iV = 4 for p = 3 , 5 and 7 with boundary parameters 
(3 + = 0.695 ,(3- = 0.774. We remark that these solutions reproduce the s = \ case given in 
[1H1I17]. The T — Q relations (13 . 1 0[) and (13.111) are also numerically verified for s = | case for 
some selected values of N , p and boundary parameters. We acknowledge that these analysis 
provide some numerical support for the completeness of the Bethe-ansatz-type equations 
derived, ( I3.20p and ( I3.2ip . and not a complete rigorous proof. We have also verified that the 
energies given in Tables 1 and 2 coincide with those obtained from direct diagonalization of 

5 Discussion 

By using a method that relies on certain functional relations that the "fundamental" transfer 
matrices, t^ ,s \n), obey at roots of unity and the truncation of fusion hierarchy, we set up a 
generalized form of the T — Q relation, (13.101) and (13. lip , for the open spin-s XXZ quantum 
spin chain with nondiagonal boundary terms. From these relations, we have determined 
Bethe-ansatz-type solutions of the model, (I3.20p and (I3.2ip . These solutions hold only for 
V — f ,f The solutions found here hold for arbitrary values of boundary parameters (at 
most two). These solutions have been checked for chains of length up to N = 4 for p = 3 , 5 
and 7 with boundary parameters (3 + = 0.695 , /?_ = 0.774. We verified that they indeed 
produce all the (2s + 1)^ transfer matrix eigenvalues for s = ~ , 1 and |. Moreover, we also 
presented numerical evidence for the completeness of the Bethe-ansatz-type solutions found 
(using s = 1 as examples) in Tables 1 and 2. The numerical support for the completeness 
of the solutions presented here (using s = 1 case as examples) together with the results 
presented for spin-1/2 case in 1 16. 47] and the fusion hierarchy ( I2.13P which is used in the 
construction of higher spin-s transfer matrices could perhaps possibly enable one to develop 
a more formal rigorous proof for the completeness of the solutions found here. It would be 
interesting to pursue this in the future. 

In addition, a number of problems remain that are worth investigating. Perhaps one 
could carry out a more thorough treatment and analysis of the functional equation to yield 
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the exact form of the Qi(u) functions, thus avoiding the need for an ansatz such as (I3.14p . 
Another interesting problem is to see the relation of s — 1 case to the supersymmetric 
sine-Gordon (SSG) model, along the lines of [58] and [59], but now for spin-1 chain with 
nondiagonal boundary terms described by the generalized T — Q relations instead of the 
conventional T — Q relation. One could also try to generalize the solutions presented in [H] 
for the spin-1/2 case, where all six boundary parameters are completely arbitrary, to any 
spin s, and analyze the s = 1 case for this general solution in relation to the SSG model. 
In this regard, one can study the continuum limit of their Nonlinear Integral Equations 
(NLIEs), thus investigating the infrared (IR) and ultraviolet (UV) limits of the NLIEs. One 
could also investigate the boundary bound states of SSG models corresponding to all these 
cases such as reported recently in [85J. We hope to be able to address some of these issues 
in future publications. 
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E 


Bethe roots, {uf } } 


-5.6483 


0.426847 + 2.19193 i, 0.719676 + 1.1781 i, 0.109151 i, 
0.426847 + 0.164266 i 


-4.67715 


0.106242 + 2.28424 i, 0.379199 + 1.1781 i, 1.05101 + 1.1781 i, 
0.106242 + 0.071957 i 


-2.75841 


0.387014 + 2.748893 i, 1.277532 i, 0.932369 + 1.1781 i, 

0.0609966 i 


-1.98286 


0.185547 + 2.748893 i, 1.701637 i, 0.915819 + 1.1781 i, 

0.138044 i 


-1.54571 


0.171807 + 3.046499 i, 0.171807 + 2.451287 i, 1.566925 i, 
0.916569 + 1.1781 i 


-0.489791 


0.781754 + 1.921787i, 1.599981 i, 0.0312436 i, 
0.781754 + 0.434407 i 


-0.392189 


3.109568 i, 0.779636 + 1.920991 i, 1.554992 i, 
0.779636 + 0.435203 i 


0.572634 


0.810472 i, 0.624212 + 1.1781 i, 0.010646 i, 
1.227343 + 1.1781 i 


0.808501 


3.130312 i, 0.791507 i, 0.618753 + 1.1781 i, 
1.221033 + 1.1781 i 



Table 1: The 9 energies and corresponding Bethe roots for N — 2,s — l,p — S,rj — 
in/ '4 , a_ = , /3_ = 0.767 , 9_ = 0.573 , a+ = , = 0.598 , 9+ = 0.573 
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E 


Bethe roots, {v£ } } 


-6.07709 


0.0471453 + 3.1415 i, 0.0471453 + 2.61809 i, 1.74867 i, 0.74532 + 1.309 i, 

0.48742 i 


-4.65604 


2.65564 i, 0.107433 + 2.35618 i, 0.321204 + 1.309 i, 0.557414 i, 
0.107433 + 0.261819 i 


-4.3506 


0.00657235 + 3.07819 i, 0.00657235 + 2.6814 i, 2.07693 i, 0.12098 + 1.93837 i, 

0.12098 + 0.679624 i 


-2.55991 


0.272597 + 3.13706 i, 0.272597 + 2.62253 i, 2.13098 i, 0.672718 + 1.309 i, 

0.862768 i 


-1.63092 


0.326829 + 2.87979 i, 0.308315 + 2.35663 i, 2.13093 i, 0.890835 i, 

0.308315 + 0.261367 i 


0.0925845 


0.248529 + 2.87979 i, 1.76311 i, 0.373083 + 1.309 i, 1.20497 + 1.309 i, 

0.487 i 


0.0971716 


0.548694 + 2.59187 i, 2.13099 i, 0.518481 + 1.309 i, 0.856853 i, 
0.548694 + 0.0261235 i 


1.6757 


0.70468 + 2.87979 i, 0.338436 + 1.309 i, 0.854426 i, 1.08306 + 1.309 i, 

0.487 i 


2.99332 


1.7639 i, 0.273003 + 1.309 i, 0.720682 + 1.309 i, 0.487 i, 
1.54847 + 1.309 i 



Table 2: The 9 energies and corresponding Bethe roots for N — 2,s — l,p — 5,rj — 
in/6 , a_ = 0.854* , /5_ = , 9_ = 0.482 , a+ = 0A87i ,[3 + = 0,9 + = 0.482 
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